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The evolution of the r-mode instabihty is hkely to be accompanied by secular kinematic effects which will 
produce differential rotation with large scale drifts of fluid elements, mostly in the azimuthal direction. As 
(— >^ I first discussed in [hh, the interaction of these secular velocity fields with a pre-existing neutron star magnetic 

field could result in the generation of intense and large scale toroidal fields. Following their derivation in the 
companion paper B, we here discuss the numerical solution of the evolution equations for the magnetic field. 
The values of the magnetic fields obtained in this way are used to estimate the conditions under which the r-mode 
instability might be prevented or suppressed. We also assess the impact of the generation of large magnetic fields 
on the gravitational wave detectability of r-mode unstable neutron stars. Our results indicate that the signal to 
noise ratio in the detection of gravitational waves from the r-mode instability might be considerably decreased 
if the latter develops in neutron stars with initial magnetic fields larger than 10^" G. 

PACS numbers: 04.70.Bw, 04.25.Dm, 04.25.Nx, 04.30.Nk 
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f — . I. INTRODUCTION 

o. 

-«,^ Hot, newly born and fast spinning neutron stars are among the best candidates for triggering an r-mode instability. These 

^^ candidates, however, are also likely to be threaded with intense magnetic fields in the range 10^^ — 10^'^ G, that are usually 
O neglected when discussing the onset and growth of the r-mode instability. As in our introductory paper [jl|], we here focus our 
^Y attention on the presence of these magnetic fields. In particular, we develop in more detail the evolutionary scenario outlined 
5-H in [|[] (hereafter paper I), where we have argued that the kinematical properties of the r-mode oscillations will produce secular 
Pf velocity fields that couple to any pre-existing magnetic field and result in a net amplification of the latter. 

J> Using the set of induction equations derived in paper I, we here present the results of their numerical integration and show that 

the exponential growth of the mode amplitude is also accompanied by an exponential growth of a toroidal magnetic field. By the 
time that the instability has reached saturation, the newly generated toroidal field might become comparable with or larger than 
j^ the seed poloidal magnetic field if the saturation amplitude is of order unity. The generation of magnetic field is accompanied 
by a conversion of the energy of the mode into magnetic energy. By computing the growth rate of the magnetic field, we can 
calculate the rate of mode energy loss to magnetic energy and estimate the strength of the magnetic field which would either 
prevent the onset of the instability or suppress its evolution. In our calculations we consider separately the case in which the 
stellar core is composed of a normal neutron fluid from the case in which the core is a Type II superconductor. 

The generation of magnetic field also decreases the amount of energy which can be radiated to infinity in the form of gravita- 
tional waves and we here discuss how the expectations for detectability of r-mode unstable neutron stars through gravitational 
radiation should be modified in terms of the initial magnetic field in the star and of the mode amplitude at saturation. 

The following is a brief outline of the contents of the paper. In Section n we present our model for the evolution of the 
r-mode instability in a magnetic neutron star. In particular, we first discuss the background evolution of the star's angular 
velocity and of the mode amplitude, and then present our calculations for the evolution of the magnetic fields. Next, we consider 
in Section |n|the impact of the newly generated magnetic fields on the evolution of the instability, estimating the critical magnetic 
fields for the prevention of the instability or its suppression. Section W^ is devoted to a revision of the standard picture for the 
gravitational wave detectability of the r-mode instability for a magnetic neutron star. Finally, conclusions and perspectives for 
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future developments are presented in Section VII 



II. EVOLUTION OF THE i?-MODE INSTABILITY IN A MAGNETIC NEUTRON STAR 

The description of the fully developed character of the instability is a non trivial task and, at present, several important aspects 
of the instability still need clarification. A self-consistent calculation of the problem requires fully nonlinear magnetohydro- 
dynamics as well as General Relativity. This is beyond the scope of this paper However, some important conclusions can be 
drawn by considering the 7--mode instability and the generation of magnetic fields as distinct processes, evolving independently. 
This is equivalent to assuming that the nature of the r-mode oscillations is unaffected by the generation of large magnetic fields 
and can, at any time, be described by the perturbative expressions [cf. Eq. (1) in paper I]. Within this framework, we can then 
specify a model for the evolution of the r-mode instability and couple it with the evolution of the magnetic field generated by 
the secular effects produced in r-mode oscillations, and discussed in detail in paper I. 



A. Numerical model of r-mode evolution 



In order to provide a direct comparison with results presented in the literature we have made use of the phenomenological 
model for the evolution of the instability presented by Owen et al. [0] in the case of a rotating neutron star with zero magnetic 
field. In this model, the evolution of the r-mode instability is assumed to have three phases. The initial phase is the one during 
which any infinitesimal axial perturbation will grow exponentially over the timescale t^^, set by the gravitational radiation- 
reaction. For an ^ = m = 2 mode|^and a neutron star initially rotating at the "break-up" limit fJ^, at a temperature ~ 10^" K 
(and for a number of different equations of state), this timescale has been estimated to be of the order of a few tens of seconds. 
The exponential growth is followed by an intermediate phase during which the amplitude of the mode is expected to reach its 
saturation value as^t, and the star is progressively spun down as a result of the angular momentum loss via gravitational waves. 
This phase has been estimated to be of the order of about one year if conventional cooling rates and viscosity estimates are used. 
The final phase of the evolution is the one in which viscous dissipative effects dominate the radiation-reaction forces and the r 
modes start to be damped out. 

Phenomenological expressions for the time evolution of the mode amplitude a and of the star's angular velocity il can be 
derived by requiring that the loss rates of the energy in the mode and of the star's total angular momentum are the same as those 
recorded at infinity in the form of gravitational waves. The ordinary differential equations governing this evolution have been 
first derived in [|[| and can be synthesized as[| 
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where r^ is the viscous timescales comprising both bulk and shear viscosity (see [^] for a definition) and Q is a constant [see 
Eq. (25) of paper I for a definition]. 

The numerical solutions to equations (|l|) and (g) are shown in Fig. [l|, where we have plotted the time evolution of the mode 
amplitude (dotted lines) and of the star's angular velocity (solid lines), normalized to the saturation value and the break-up limit 
respectively. The main panel of Fig. |l] shows the initial phases of the amplitude growth and angular velocity decay and is, for 
this reason, shown on a linear temporal scale. Different solid curves refer to different values of the saturation amplitude. Note 
that the decrease in the star's angular velocity and that its value after one year depends sensitively on the value used for agat and 
becomes very small only for agat = 1- This is apparent from the inset which shows the decay of the star's angular velocity on a 
(longer) logarithmic time scale. 



^Hereafter we will only consider modes with £ = m. 

■^Note that our expression for dVl/dt after saturation differs from the one given in [B], where it was assumed that Oa 
always 0(1), independently of the actual value of Qsat. 
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FIG. 1. Time evolution of the mode amplitude (dotted lines) and of the star's angular velocity (solid lines) normalized to the saturation 
value Qsat and the break-up limit fi^ , respectively. Different curves refer to different values of the saturation amplitude. The small inset shows 
the decay of the star's angular velocity on a logarithmic time scale. 



Next, we consider the effect of an increasing mode amplitude on the kinematics of r-mode oscillations. For doing this we 
compute the equations of motion allowing the mode amplitude to vary as in Eq. (lit) in the case before saturation. Fig. shows 
the effect of a variable mode amplitude on the trajectories of two fiducial fluid elements following r-mode oscillations with 
£ — 2. We concentrate on two fluid elements which have the same initial longitude (J)q — and latitude ^ 50 degrees, with the 
mode having an initial amplitude ao — 0.02, which then saturates at asat = 0.1. The initial positions are indicated with stars, 
while the positions at the end of each oscillation are indicated with filled squares. Note that in order to have a growth time which 
is larger but comparable with the oscillation period, we have artificially set O = bfl^ and the trajectories shown in Fig. g should 
therefore be considered only as indicative. 

When compared with the coiTesponding trajectories shown in Fig. n^ of paper I, they show that the effect of an exponentially 
increasing mode amplitude is that of exponentially increasing the azimuthal and polar drifts between oscillations. As discussed 
in paper I when deriving the expression for the azimuthal drift [cf. Eq. (6) of paper I], one should expect secular effects to emerge 
every time that the velocity field at the position of the fluid element varies as it moves. This can happen either because the fluid 
velocity has a dependence on the position of the fluid element (as in the case of constant amplitude oscillations) or because 
the amplitude of the fluid velocity is changing (as during the growth of the instability). As a result, a mode amplitude which 
is exponentially growing in time will produce a secular motion, independently of whether the modes have reached nonlinear 
amplitudes. 
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FIG. 2. Projected trajectories 9 (t) sin 6 (t) cos <j)(t), and <j)(t) sin 9{t) cos (l>{t) of two fiducial fluid elements as seen in the corotating 
frame. While the fluid elements are allowed to oscillate, the amplitude of the mode increases from an initial value ao — 0.02 to a saturated one 
cvsat = 0.1. The initial positions are indicated with stars, while the positions at the end of each period are indicated with filled squares. Note 
that we have set SI — Sfi^ ; this is to have a growth time which is larger but comparable with the oscillation period and to show schematically 
the effects of a increasing mode amplitude on the secular motions. 

In the following Section, where we discuss the numerical evolution of the magnetic field as a result of the secular velocity 
field, it will be possible to distinguish the contribution to the field generation coming from the exponential growth phase of the 
mode from the one resulting from the secular drift at mode saturation. 



B. Evolution of the stellar magnetic field 



We here present results from the solution of the induction equations as obtained with both the Lagrangian method [cf. Eq. 
(14) of paper I] and with the orbit-average method [cf. Eq. (16) of paper I] discussed in paper I. The computations have been 
performed for a number of different values of the saturation amplitudes, for modes with i = 2,3, 4, and with an initial dipolar 
magnetic field [cf. Eq. (8) of paper I]. 

From a numerical point of view, the Lagrangian and orbit-average approaches have a common strategy. Both of them, in 
fact, solve the equations of motion [cf. eqs. (3) of paper I] for a number of fiducial fluid elements suitably distributed on an 
isobaric surface of the star and calculated for several different mode numbers. In the Lagrangian approach, in particular, the new 
positions after each oscillation are used to reconstruct the strain tensor Sjk — dx^ {t)/dx^ (to) from which the new magnetic field 
components are computed. Conversely, in the numerically more expensive Eulerian approach, the secular polar and azimuthal 
velocities v^ , v'^ are computed after each oscillation period and their values transformed onto a finite difference grid corotating 
with the star, over which the induction equations are then solved. In this latter case, a second-order Lax-Wendroff evolution 
scheme is used to minimize artificial dissipative effects and accurately estimate the field growth rate. 
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FIG. 3. Growth of the logarithm of {AB'^) scaled to the initial poloidal field. The left panel refers to an ^ = 2 mode and a number of 
different values of Osat. The right panel shows the growth rate for different £ modes and two different saturation amplitudes. The underlying 
star has an initial angular velocity Slo equal to the break-up limit fl^ , and (-Bo) = 10^ G. Note that in plotting the evolution of different mode 
numbers in the right panel we have assumed they all have the same growth time. This allows us to present them on a single plot but one should 
bear in mind the higher-order modes have considerably larger growth times ||3|. 



Indicating by (AS*) the i-th component of the volume average secular magnetic field, we show in Fig. Bthe time evolution 
of the averaged toroidal magnetic field normalized to the average initial magnetic field. The left panel of Fig. ^ in particular, 
shows the time evolution of the average toroidal field as computed with the Lagrangian approach (continuous lines) and with 
the orbit average approach (dotted lines). Different lines refer to different values of the saturation amplitude ckgat and have been 
computed for an i — 2 mode and a "fiducial" neutron star with mass M — 1.4M0, radius R — 12.5 km, and initial angular 
velocity flo equal to the "break-up limit" fig = {2/3)y/T:Gp, with p — 3Af/(47r_R'^). The value of the initial toroidal field 
served only to set a lower limit to the vertical axis and was here set to be 10^^ (Bq). Each line in the left panel shows an initial 
stage during which the toroidal field grows exponentially and a subsequent one during which the field has a growth which is a 
power-law of time. As mentioned in the previous Section, the two different regimes in the field growth reflect the fact that there 
are two regimes in which secular kinematic effects become apparent. 

The first regime is related to the exponential growth of the mode amplitude and produces most of the new magnetic field. 
During this phase, the secular effects are extremely large [i.e. ©(a^ exp(t))] and the toroidal magnetic field is either produced 
or amplified by the wrapping of the poloidal magnetic field produced by the (mostly) toroidal secular velocity field. The 
amplification is so large that the toroidal magnetic field soon becomes comparable to and then larger than the seed poloidal field. 
Larger values of the saturation amplitude produce proportionally larger amplifications and these can be so dramatic in the case 
of ttsat — 1, that the newly generated magnetic field at mode saturation has become about four orders of magnitude larger than 
any pre-existing magnetic field. Note that a magnetic configuration with strong toroidal magnetic fields is generally unstable 
and hoop magnetic stresses in the radial direction will tend to move magnetic field lines towards the poles [^]. Although the 
instability of this configuration might lead to a number of interesting astrophysical phenomena (see, for example, [Q] where the 
buoyancy instability of this configuration is used for a 7-ray burst model), we will here assume that gas pressure gradients are 
always be able to counteract these hoop stresses. 

The second regime is instead related to the saturated stage of the transition, during which the mode amplitude is constant in 
time. During this phase, the secular effects are much smaller [i.e. 0{a^^t.t)] and so is the field growth rate, which is a power-law 
of time. 

The left panel of Fig. also gives a comparison between the Lagrangian and orbit-average methods. Despite the differences 
between the two approaches, the agreement is extremely good and deviations between the two appear only for the late-time, 
ttsat = 1 growth rate, and at the beginning, when the growth rate is extremely small. This difference appears because of the 



difficulty of calculating the particle trajectories accurately for such large values of the mode amplitude and the inadequacy of the 
Lagrangian method in this regime. 
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FIG. 4. Time evolution of (AS"") (left panel) and of (AS*) (right panel) scaled to the initial field strength. The calculation has been done 
with the orbit-average method for an i" = 2 mode and a saturation amplitude Osat ~ W~^. 



The right panel of Fig. B^ shows the same quantity as the left one, but for different mode numbers and two different saturation 
amplitudes. Note that, in general, r modes with different mode numbers saturate on different timescales, with the high-order 
modes having considerably longer growth times. However, in order to present them all in the same plot, we have assumed in 
the right panel of Fig. ^that they all have the same timescale. Interestingly, higher mode numbers would, in principle, be more 
efficient in generating a toroidal magnetic field, and thus in converting the energy of the mode into magnetic energy, but modes 
with £ > 2 are also very unlikely to become unstable because of their longer growth times. 

In Fig. we show the corresponding evolution of the poloidal magnetic field components for an £ = 2 mode and a saturation 
amplitude agat ~ 10^^ as obtained through the orbit average method. Note that neither of the two poloidal magnetic field 
components shows a significant secular growth and that they both oscillate around the initial values. This is in agreement with 
Bragisnky's dynamo model, whose velocity field has many analogies with the secular velocity field set up by the r modes, and 
which predicts that the poloidal magnetic field should not be subject to a dynamo action and have a zero time average (see [||] 
for a complete discussion of Braginsky's dynamo). Fig. H basically states that a large toroidal magnetic field will be the main 
outcome of the secular kinematic effects produced during the r-mode instability. 

The numerical integration of the equations of motion and the calculation of the magnetic fields produced was usually stopped 
after 2000 s, corresponding to ^ 10^ r-mode oscillations. At that stage and for our standard neutron star, the instability 
has already reached saturation and the toroidal field growth has entered a power-law growth phase which can be described 
analytically (cf. Fig. ||). 

Fig. S synthesizes the results presented in the previous Figures and shows the time evolution of the total magnetic field (AB) 
scaled to the initial magnetic field and over a timescale of one year The left panel, in particular, refers to an £ = 2 mode, 
an initial angular velocity Oq = fl^ and three different saturation amplitudes. The panel also presents a comparison between 
the fully numerical Lagrangian solution and the analytical one based on equations (26) and (27) of paper I. The agreement is 
remarkably good even for agat — 1 when we expect the expansion procedures in powers of a to be rather poor. 
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FIG. 5. The left panel shows the long-time evolution of the logarithm of total magnetic field (AB) scaled to the initial field for an £ = 2 
mode and different saturation amplitudes. The star is set to have Qo ~ Sl^ . The right panel is the same as the left one but shows the different 
evolutions resulting from values of the initial angular velocity Qo = jQ^ , fio ~ ^^b > ^o = |^b , ^o = ^b ■ 



The right panel of Fig. g is the same as the left one but also shows the effect of a smaller initial angular velocity (i.e. 
fio — j^B J ^0 = 5^B 7 ^0 = j^B I ^0 — ^b)- Clearly, for stars rotating at slower angular velocities, the instability and the 
field growth set in at progressively later times. However, for any given saturation amplitude, the final magnetic field produced 
is rather insensitive of the rotation rate. This is because the magnetic field is generated mostly during the phase of exponential 
growth of the mode and as long as the star rotates sufficiently fast for the mode to grow unstable, the end result will remain 
roughly unchanged. 



III. IMPACT OF MAGNETIC FIELDS 



In paper I, as well as in the previous Sections of this paper, we have described how r-mode oscillations can produce secular 
drift velocity fields which, in a highly conducting plasma such as neutron star matter, can interact with pre-existing magnetic 
fields producing new ones. We have also discussed the set of equations that can be solved to calculate the magnitude of the newly 
generated magnetic fields and have solved them numerically to compute the values of the different magnetic field components 
produced for a number of different mode amplitudes and stellar angular velocities. We have not yet discussed the impact of 
these fields on the existence or growth of the r-mode oscillations. As mentioned in Section ||, a fully self-consistent treatment 
is beyond the scope of this paper However, important considerations can be made if we treat the r-mode instability and 
the generation of magnetic fields as evolving independently and with the magnetic field not "back-reacting" on the kinematic 
features of the r-mode oscillations. Under this simplifying hypothesis, it is possible to calculate the strength of the magnetic 
field necessary to prevent the amplification of the first r-mode oscillation, or damp the instability when this is free to develop. 
Both of these aspects are discussed in the following Sections IV and Q. 



IV. DISTORTION OF R MODES BY THE STELLAR MAGNETIC FIELD 



As the magnetic field strength increases, the magnetic back-reaction forces generated by the r-mode motions of the fluid 
become larger and larger As the magnetic field strength increases, the normal modes of the star coiTesponding to the zero- 
magnetic-field r modes differ in character more and more from the zero-magnetic-field r modes, taking on appreciable hydro- 
magnetic and hydromagnetic-inertial properties. Eventually, the velocity perturbations corresponding to these normal modes are 



significantly different from the zero-magnetic-field r modes to the point of preventing gravitational radiation from amplifying 
these distorted r-mode oscillations. 

Since the normal modes of the star are global in character, the effects of a magnetic field cannot be analyzed fully by a local 
analysis, such as can be done by making a local comparison of forces. In contrast, a mode-energy approach takes into account 
the global properties of the modes and magnetic field and this is what we will consider in the following. During an oscillation the 
energy in the magnetic field will rise and fall. If there is not enough energy in the mode to supply the maximum magnetic energy 
increase required to complete an oscillation, a "full" r-mode oscillation will not occur (some other, smaller scale pulsations 
might still occur). Stated differently, if the magnetic field exceeds a critical value, -Bciit.p, the magnetic stress that builds up 
during an oscillation will be so large that it will halt the fluid motion involved in the oscillation, i.e. the fluid momentum density 
will be brought to zero. The condition that defines Sciit.p is therefore 

SE,, ^E, (3) 

where E the energy in the mode as measured in the corotating frame of the equilibrium star 

E=- ^ /"" p r^dr ~ 8.2 X IQ-"" aHm'' R^" , (4) 







and 5E^ is the change in magnetic energy density 

5E^^ = ^ f SB^ d^^ , (5) 

with SB^ — (SBP)'^ + [SB'^Y. Note that the fractional change in the mode energy produced by gravitational wave emission 
during a single oscillation period is Att/lo\t^^\. This is always ^ 1 and therefore can be neglected in this comparison. 

The expression for 5E^ varies according to whether the neutron star core is composed of a normal neutron fluid (N) or is 
superconducting (SC). More particularly, we denote as &E^ and -B^jt „ the magnetic energy and the critical, volume averaged, 

SC SC 

magnetic field for a normal core. Similarly, we denote as oE^ and iJ^rit p ^^ corresponding quantities for a superconducting 
core. The two cases will be considered separately below. 

A. When the stellar core is normal 

The induction equation can be integrated in time to estimate the variations in the magnetic field produced by the perturbation 
velocity field as 

5B^ ~B^ f 9{t')dt' ~B^ f ^^dt' ~ SBP , (6) 

SB^ ~ bI I 4>{t')dt' ~ B'a I ^^dt' , (7) 

where we have set ^ ~ ~ 6v/r. Because of the periodic variation of the magnetic energy during an oscillation, expression (Bl) 
should be averaged over half of an oscillation period, so that the magnetic energy produced is 
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In deriving (ph we have not included contributions to SE^ from changes in the magnetic field outside the star which are heavily 
suppressed by the high electric conductivity of the crust. Also, we have written the time average of the velocity perturbation as 

;"-«'>-'== (I) I- 

Note that both energies (0) and (ph depend quadraticaUy on the velocity perturbation Sv and therefore on a. As a result, the 
critical fields are not not sensitive to the amplitude development of the mode. The volume integral can be easily performed after 
assuming an initial dipolar magnetic field: Bq — Bd{R/r)'^'9{0), where ^(0) = 4cos^ 9 + sin^ 6. Indicating with |Y^(0, 0)| 
the modulus of the vector spherical harmonic, the magnetic energy produced over half of a period is then 
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It is worth noticing that we can here neglect the distortion of the magnetic field near the crust-core boundary during an r- 
mode oscillation. The reason for this is that, on the timescales relevant for the r modes and regardless of whether the core is 
superconducting, the magnetic field Unes (and flux tubes) are pinned in place at the crust-core boundary by the high conductivity 
of the crust. 

Using (|iy) with p = 0.5 and ([ll|), we can now write the condition for r-mode prevention (g|) in terms of a critical field given 
by 



^crit.p 



2.5 X 10 



16 



B J 



G 



(12) 



where Mia = M/IAMq, R12.5 = i?/12.5km, and n*g = 5.62 x 10^ rad s^^. 



B. When the stellar core is superconducting 

When, as a result of the star's cooling, the core temperature falls below a critical value, Tc ~ 3 x 10^ K, the core is expected 
to become a Type II superconductor [jl^. In such a superconducting core, the magnetic field is confined to flux tubes and the 
total energy of the flux tubes can be computed in terms of the total number N of flux tubes in the core and of the mean total 
energy of a flux tube 



with £f being the average length of a flux tube, Hf the magnetic field inside it and 0o = hc/{2e) 
quantum. The total magnetic energy of the flux tubes of the unperturbed star is then 
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where 0tot ~ ttR^B^ is the total magnetic flux. The r-mode oscillations will shear and stretch the flux tubes, increasing their 
length by S£, and thus changing their energy content by an amount 5E^^^ = {(ptot /8Tr)Hf Si f over half a period. Estimating the 
volume-averaged variation in the flux tube length as 
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we then obtain 
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(15) 



(16) 



where we have taken Hj w H^-^ , the lower critical field [g|. 

Note that expression (|l6|) is only valid when the spacing between the flux tubes is much larger than the London penetration 
length, i.e. for (B) ^ H^-^ . However, as (B) rises above H^-^ and tends to the upper critical field H^^ w 10^^ G, Hf also rises 
and tends to H^,^ , with the flux tubes crowding one another When (B) ~ H^^ , the normal cores of the flux tubes touch, the 
superconducting state is destroyed [^ and the variation in the magnetic energy is then again given by expression (|l0|). Combining 
(llq) with (ttOb, we can now write the condition for r-mode prevention in terms of the critical field for a superconducting core 
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As for the critical field for a normal core, expression (111) is insensitive to the values of the mode saturation amplitude, underlining 
that if Bq > Scrit.p even very small amplitude fluctuations of r-mode character cannot occur. On the other hand, the critical 



fields are sensitive to the star's spin rate; in particular, if fi is w J7|j, the magnetic field needed to prevent oscillation of the 
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crit,p 



B. 



N 



^.jij , whether or not the core is superconducting. If instead il ^ ilg, then 

10^^ (fi/0.01 n*j^)^H^^QMiARi2.5 G' where Hcim = -ffci/10^^ G. 
The critical fields (n2b and (HTh are several orders of magnitude larger than the typical magnetic fields of young neutron stars 
that observations indicate to be in the range 10^^ — 10^"^ G. This suggests that the distortion of the first r-mode oscillations 
by magnetic fields is rather unlikely and that the r-mode instability in newly born neutron stars (or any compact object with 
analogous properties) is, in general, free to develop. The onset of the instability, however, will also lead to the production of 
very intense magnetic fields which could ultimately lead to the damping of the oscillations, as discussed below. 



V. DAMPING OF R MODES BY THE STELLAR MAGNETIC FIELD 

The damping of oscillations by the growing magnetic field reflects the question of whether the energy required to continue 
the oscillation for the next time step At exceeds the energy that can be pumped into the mode in the time At by the emission of 
gravitational radiation. 



A. When the stellar core is normal 



A critical point in the balance between the energy spent in producing magnetic field and the energy in the mode provided by 
the emission of gravitational waves is reached when the two rates are equal 



dt 




(18) 



where 




327rG 


4 
i3(5)!!j 



i?2 



pr dr 



(19) 



is the rate of energy transferred to an ^ = 2 mode ^]. After the equality (|l8|) is reached, dE^/dt > (dE/dt)^^. The reason 
for this is that the toroidal magnetic field, and hence dE^^/dt, will continue to grow whereas (dE/dt)^^ remains the same or 
decreases as gravitational radiation causes the star to spin down. The only source of energy to feed E^^ is thus the energy E of 
the mode which damps on a timescale 



E 



" (dE^Jdt) 
A measure of the relative importance of the two energy rates can be obtained through considering their ratio 



9{t) 



dE^Jdt 



[dE/dt)^^ 
where, in the case of normal core, the rate of production of magnetic energy takes the form 



dt 



4(1 -p) 



Qnp 



BlR^A'a^{t)n{t) / a^(t')n{t')dt' 





(20) 



(21) 



(22) 



with now A' = /p d(t)J^iK2{9))^\Y§2i^, </')l*(6') sm9d9, which is again 0(1). 

Fig. ^ shows the numerical computation of the ratio (gTj) for a "fiducial" neutron star initially rotating at the break-up limit 
fi^. The rate of energy loss to gravitational radiation has been computed for an £ = 2 mode and a poly tropic equation of state 
with r = 2. The left panel shows the time evolution of g(t) for saturation amplitudes asat =0.1 and four different values of 
the initial magnetic field (Bq) which are indicated with different types of thick lines. The right panel is the same as the left 
one except for having Ofgat — 1- For each curve in Fig. ^ it is straightforward to distinguish the initial phase of exponential 
growth of the magnetic field (and of the magnetic energy production rate) from the subsequent phase of power-law growth. It 
also appears evident that, for a given saturation amplitude, different values of (Bq) have only the effect of moving the curves in 
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the vertical direction: the smaller the initial field, the longer it will take to reach a sufficiently strong magnetic field necessary 
for suppression. 

As discussed in paper I, we expect the drift of fluid elements given by the velocity field Vd [cf. Eq. (7) of paper I] to be 
qualitatively correct to 0{a^). However, within our approach we cannot exclude the possibility that the actual drift velocity 
is smaller (or larger) than the estimated one. To assess the variations brought about by a smaller secular velocity field, we use 
the thin lines in Fig. ^ that are near to the thick lines of the same type to show the magnetic field evolution which has been 
calculated with the secular velocity field being one tenth of that computed in paper I. Because of the exponential growth in the 
mode amplitude and in the toroidal magnetic field, the overall differences produced are rather small. 
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FIG. 6. Time evolution of the ratio g{t). The calculation refers to our "fiducial" neutron star initially rotating at the break-up limit Q.%. 
The left and right panels show the evolution for different values of the saturation amplitude. Lines of different type refer to different values of 
the initial magnetic field (Bo)- Nearby thick and thin lines of the same type are used to compare the magnetic field evolution produced by the 
secular velocity field calculated in paper I [cf. Eq. (7) of paper I] and a velocity field which is just a tenth of it. 



When read in terms of the times at which suppression of the oscillations begins. Fig. indicates that for an initial magnetic 
field (i?o) = 10^" — 10^^ G, the instability starts to be suppressed after a time between 2 days and 15 minutes for a saturation 
amplitude asat ~ 0.1 (see also Table I). When this happens the total magnetic fields produced are in the range lO^"* — 10^^ G. 
The condition for oscillation suppression (jlsb can also be used to calculate ficrit: the critical spin rate above which r modes are 
excited and below which they are damped. Values of the critical angular velocity for a normal core normalized to the break-up 
angular velocity Q,*g are indicated in Table I. The Table reports the critical spin rates as a function of the initial magnetic field 
strengths for a number of different saturation amplitudes. For each pair of these parameters, we also indicate in brackets the 
time (in seconds) icrit at which the condition wm is reached. As we will discuss in Section pU the knowledge of the critical 
angular velocity is in fact very important as it provides the final (and maximum) gravitational wave frequency at which the last 
gravitational radiation is emitted. Together with the initial star's spin frequency, the maximum frequency determines the domain 
over which a signal-to-noise ratio can be calculated [cf. Eq. (pq)]. 

It is worth underlining that even when the r-mode instability is rapidly suppressed by generation of intense magnetic fields, 
there is sufficient time for the star to loose a good fraction of its rotational energy to gravitational waves. This is because the loss- 
rate of angular momentum to gravitational waves is a steep function of the star's angular velocity [cf. Eq. (0)] and is large mostly 
soon after the amplitude saturation. Therefore, for initial magnetic fields of the orderof 10^° G, most of the astrophysical consid- 
erations about the importance of the r-mode instability in spinning down newly-born neutron stars [fTj are basically unchanged. 
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N 

n^it/^s, (tcritins) 


(Bo) = 10* G 


(Bo) = 10"' G 


(Bo) = 10" G 


(Bo) = 10" G 


Qsat = 1.0 

asat = 10-1 
asat = 10-2 


0.318 (2.90 X 10^) 


0.467 (3.34 X 10=*) 
0.535 (1.37 X 10^) 


0.948 (5.24 X 10^) 
0.975 (9.75 X 10^) 
0.974 (5.62 X W) 


0.999 (3.52 X 10^) 
0.999 (3.23 X 10^) 
0.999 (3.23 X 10^) 



TABLE I. Normalized critical angular velocities [i.e. final values of the star's angular velocity when g{t) = 1, normalized to the initial 
break-up angular velocity Vlg] for different values of the magnetic field (Bo) and saturation amplitude Qsat. Values not reported refer to times 
longer than one year, for which the instability is suppressed by viscosity rather than by the production of magnetic field. The numbers in 
brackets indicate the time (in seconds) after which the condition for suppression is reached. 



B. When the stellar core Is superconducting 

In the case of a superconducting core, the calculation of the damping of the r mode oscillations is slightly more complicated 
(and uncertain) because additional considerations need to be made about the time variations of the total magnetic flux 0tot and 
about the magnetic field components that contribute mostly to the flux. Once B'^ exceeds B'', the magnetic energy rate can be 
expressed as 



dE^ (t) ., ^2(9) 



H^^BPR^a^{t)n{t) 



dt 60 

Hence, for a superconducting core, the actual critical spin rate is likely to be 



(23) 



^. 



crit 



O.Tffi 



BP 



^ I 1012G 



1/7 



H 



f 



1016G 



1/7 



7,r-2/7p-3/7 



(24) 



which is again independent of the mode amplitude. 



VI. GRAVITATIONAL WAVE DETECTABILITY 



The relevance of the gravitational radiation emitted by r-modes oscillations has been investigated by a number of authors who 
have considered them both as single sources [| 1 2 , 13[] and as a cosmological population producing a stochastic background Jl4| ] . 
Their estimates indicate that sufficiently close sources might be detectable by LIGO and VIRGO when they reach their "en- 
hanced" level of sensitivity or that two nearby interferometers with "advanced" sensitivities might be used for the detection. 

Using the results obtained in the previous Section and bearing in mind the simplifying assumptions under which they have 
been derived, we now re-evaluate the possibility of detecting gravitational waves from r modes in magnetic neutron stars. All 
of the results discussed below refer to a "fiducial" neutron star with a normal core. For an £ = 2 mode, the rate of angular 
momentum loss due to the emission of gravitational waves can be written as [pt] 



dJc D' 2 



hi)) 



(25) 



where h^ and h^ are the strain amplitudes of the two polarizations of gravitational waves, D is the distance to the source, and 
((. . .)) denotes the average over the orientation of the source and its location on the sky. Here Jc — —SJAIR^flga^ is the 

canonical angular momentum JI^ ] (see also S) and J — Jq pr^dr/{MR'^). 
FoUowing ||^], the averaged gravitational wave amplitude is then given by 



Kt)^\lYQ^^K 



K))-- 

""" 135 



TT GM an^R^ J 



(26) 



12 



and its time evolution for a "fiducial" neutron star placed at 20 Mpc is shown in Fig. 0. The left panel, in particular, shows 
evolutionary tracks for a number of different initial values of the magnetic field and a saturation amplitude agat = 1- The tracks 
terminate when g{t) — 1, which is marked with an open symbol (different symbols refer to different initial magnetic fields). 
The right panel, on the other hand, refers to an initial magnetic field (Bq) = lO^*' G and shows the effect of different saturation 
amplitudes on the reaching of a magnetic energy production rate sufficient to suppress the instability (which occurs at the point 
marked with filled squares). Taking a saturation amplitude agat = 0.1, one sees that the emission of gravitational waves would 
start being suppressed after about 20 days. 
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FIG. 7. Time evolution of the averaged gravitational wave amplitude. The left panel is for a saturation amplitude Qsat = 1 and shows 
the different behaviour with different initial values of the magnetic field. The evolutionary tracks terminate when g{t) = 1 and different final 
points are marked with different open symbols. The right panel refers to an initial field of (-Bo) = 10^^ G and shows the influence of the 
saturation amplitude on the evolution of the wave amplitude. A filled square indicates the time at which g{t) — 1. 

A more useful quantity for discussing the detectability of the gravitational waves emitted is the characteristic amplitude he, 
defined as 



h,if)^h{t)Jf 



(27) 



where / is the gravitational wave frequency and / = cj/(27r) = 2r2/(37r) for a £ = 2 mode. 

In Fig. ^we show the frequency evolution of the characteristic amplitude for a magnetic neutron star located at 20 Mpc and 
with a saturation amplitude asat = In As in Fig. n[ different symbols mark the frequency at which g{t) — 1 for different values 
of the initial magnetic field. The frequency evolutionary track should be compared to the expected root-mean-square strain noise 
hrms curves, indicated with dashed lines Fig. and refeiTing to the LIGO "first interferometers" (LIGO I), to the "enhanced 
interferometers" (LIGO II), and to the "advanced interferometers" (LIGO III). Analytic expressions for the expected /ij-ms have 
been presented in [Bfl . 



^Because h{t) oc Qsat, but \dt/df\ oc a^^^, the characteristic amplitude is basically insensitive to Qs^ 
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FIG. 8. Frequency evolution for the characteristic amplitude for a magnetic neutron star located at 20 Mpc and with a saturation amplitude 
Osat ~ 1- As in Fig. pi different symbols mark the frequency at which g(€) = 1 for different values of the initial magnetic field. The three 
dashed lines show the expected /irms for the three stages of development of the LIGO interferometers. Note that the spike in he is the result of 
a vanishing df /dt at the onset of the instability and does not contribute significantly to the SNR. 



The characteristic amplitude also allows a direct calculation of the signal-to-noise ratio (SNR) for a single source at a distance 
D. Using matched filtering techniques, the SNR can be written as [B[| 



^1 ^2 

N 



U 



f 



(28) 



where, for a strain noise with power spectrum Sh{f), ^rms = \/fShif)- The integral in (^Sj) is computed between the minimum 
frequency /min = 251o/(3vr) and the maximum frequency /max = 2ricrit/(37r) at which the gravitational radiation is emitted, 
but is largely dominated by the amplitude-saturated part of the evolution. As a consequence, it is strongly influenced by the 
angular velocity at which the mode is suppressed ilait and, as discussed in the previous Section, on the initial strength of the 
magnetic field. 

In Table II we report the results of the calculation of signal-to-noise ratios obtained through the evolutionary models discussed 
in the previous Sections. The SNR are given for agat = 1, (o^sat = 10~^) and different values of the initial magnetic field. For a 
zero initial magnetic field and Ogat = 1, we recover the results obtained by [Bfl 
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S/N : Qsat = 1 

(a.at = 10-^) 


(Bo) = G 


{Bo) = 10» G 


(Bo) = 10^" G 


(Bo) = 10^2 Q 


(Bo) = 10" G 


LIGO-I 
LIGO-II 
LIGO-III 


1.2(0.7) 
7.6(4.1) 
10.6 (5.4) 


0.8 (— ) 
4.8(4.1) 
6.4 (5.4) 


0.4 (0.4) 
2.2(2.1) 
2.9 (2.7) 


0.1 (0.1) 
0.8 (0.6) 
1.1 (0.8) 


< 0.01 
<0.01 
<0.01 



TABLE II. Signal-to-Noise ratios for different values of the saturation amplitude asat and of the initial magnetic field (Bo). The values 
referring to a zero initial magnetic field coincide with those calculated by pi]. As for Table I, the values not reported refer to times longer than 
one year. 

Note that the presence of a relatively weak initial magnetic field (i.e. (Bo) — 10^ G) does not significantly change the SNR 
as compared to an neutron star with zero magnetic field. However, this is not the case when the initial magnetic field is taken to 
be stronger and, in particular, for an initial (Bq) — 10^^ G, the SNR is such that even LIGO-III would only be able to detect it 
if ftsat is C(l). Overall our results indicate therefore that the development of the r-mode instability in a magnetic neutron star 
could significantly reduce the expectations of detection of gravitational waves from r-mode oscillations. 

VII. CONCLUSIONS 

We have computed numerically the growth of the magnetic field generated by the secular kinematic effects emerging during the 
evolution of the r-mode instability. We have treated the magnetic field growth as a kinematic dynamo problem, thus neglecting 
any back-reaction of the magnetic field on the forces driving the velocity field. The toroidal magnetic field computed in this way 
exhibits an exponential growth which is a direct consequence of the exponential development of the r-mode instability. During 
this stage of the growth, the toroidal field can become extremely large and, in the case of a saturation amplitude of order unity, it 
can grow to be four orders of magnitude larger than the seed poloidal field. At mode saturation and thereafter, the field growth 
rate is smaller but nonzero, thus still contributing significantly to the production of very intense magnetic fields. 

Although we have neglected the back-reaction of the magnetic field on the kinematics of the r-mode oscillations, we can 
estimate the strength of the magnetic fields that would either prevent the onset of the r-mode instability or damp it when this 
is allowed to grow. Our estimates indicate that the critical initial magnetic field which would suppress even the first r-mode 
oscillation is (i3o)crit ^ 10^^ G and this is basically independent of whether the core is superconducting or not. Such a magnetic 
field is two orders of magnitude larger than those usually associated with young neutron stars and it is therefore likely that most 
newly born neutron stars will have an initial magnetic field which is sufficiently weak for the instability to be triggered. However, 
depending on the strength of the initial magnetic field and on the mode amplitude at saturation, an r-mode instability which is 
initially free to develop could be subsequently damped by the generation of magnetic fields on timescales much shorter than the 
ones set by viscous effects. As an example, with an initial magnetic field strength (Bq) — 10^" G and a normal stellar core, we 
find that the instabiUty would be suppressed in less than four hours for a saturation amplitude ckgat = 1 and in less than ten days 
forasat = 0.1. 

An important point to stress is that even in the case the r-mode instability is rapidly suppressed by the coupling of the 
mass-currents with magnetic fields, the star will have had sufficient time to loose a good fraction of its rotational energy to 
gravitational waves. This is because the loss-rate of angular momentum to gravitational waves is a steep function of the star's 
angular velocity and is therefore most efficient only in the very initial stages of the instability. As a result, for realistic values 
of the initial magnetic field, most of the astrophysical considerations about the importance of the r-mode instability in spinning 
down newly-born neutron stars remain basically valid. 

The validity of our scenario could be verified in a number of ways. As discussed in the previous Section, the larger the initial 
magnetic field, the sooner the condition for mode damping will be reached. We could therefore expect a positive correlation of 
the stellar spin rate at about one year and the initial poloidal magnetic field of the neutron star For the same reason, we could 
expect an anti-correlation between the stellar spin rate and the toroidal magnetic field at about one year. Our results also have 
a direct impact on the gravitational wave detectability of r modes from unstable neutron stars. Since the SNR for the detection 
of gravitational waves from r-mode oscillations depends sensitively on the star's spin frequency at mode suppression, and the 
latter is determined by the mode amplitude at saturation and the initial magnetic field, we have found that the SNR could be 
significantly above unity only if the initial magnetic field is below 10^'^ G. Moreover, because of the quantitative difference in 
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the mode damping in the case of a superconducting stellar core, the determination of the critical frequency could be used to infer 
a signature of the onset of proton superconductivity in the core. 

Finally, we note that some fundamental issues remain unresolved. A most outstanding one is whether the magnetic field will 
distort the r-mode velocity field in such a fine-tuned way so as to cancel the fluid drift before it suppresses the oscillation. As 
discussed in nin, considerations of the energy densities in play suggest that this is rather unlikely, but a more detailed analysis of 
this would certainly provide important information about the evolution of the r-mode instability in magnetic neutron stars. 
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